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ABSTRACT 

We report the first extrasolar planet observations from the 10.4-m Gran Telescopio 
Canarias (GTC), currently the world's largest, fully steerable, single-aperture optical 
telescope. We used the OSIRIS tunable filter imager on the GTC to acquire high- 
precision, narrow-band photometry of the transits of the giant exoplanets, TrES-2b 
and TrES-3b. We obtained near-simultaneous observations in two near-infrared (NIR) 
wavebands (790.2 and 794.4±2.0 nm) specifically chosen to avoid water vapor absorp- 
tion and skyglow so as to minimize the atmospheric effects that often limit the precision 
of ground-based photometry. Our results demonstrate a very-high photometric pre- 
cision with minimal atmospheric contamination despite relatively poor atmospheric 
conditions and some technical problems with the telescope. We find the photomet- 
ric precision for the TrES-2 observations to be 0.343 and 0.412 mmag for the 790.2 
and 794.4 nm light curves, and the precision of the TrES-3 observations was found 
to be 0.470 and 0.424 mmag for the 790.2 and 794.4 nm light curves. We also dis- 
cuss how future follow-up observations of transiting planets with this novel technique 
can contribute to the characterization of Neptune- and super-Earth-size planets to be 
discovered by space-based missions like CoRoT and Kepler, as well as measure at- 
mospheric properties of giant planets, such as the strength of atmospheric absorption 
features. 

Key words: planetary systems - stars: individual (TrES-2, TrES-3) - techniques: 
photometric 



1 INTRODUCTION 

Space-based missions like Co RoT and Kepler have the pho- 
tometric capability to detect transiting planets with radii 
not much larger t han the Earth ' s, commonly referred t o 
as "super-Earths" l|Borucki et al.ll2009l : ILlger et al.l |200ST ). 
While space-based observatories have provided the high- 
est photometric precisions, ground-based follow-up obser- 
vations play an essential role in confirming detections of 
transiting planets and characterizing the planets' orbits, in- 
teriors, composition s and atmospheres ([Torres et al.ll2008l : 
I Johnson et al.| [2009). In particular, because CoRoT and Ke- 
pler's white light observations are affected by stellar limb 
darkening (LD), there can be strong correlations between 
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the stellar LD parameters, the transit impact parameter, 
and the transit duration (which depends on the eccentricity 
and pericentre direction). High-precision photometry con- 
ducted at far-red, near-infrared (NIR) or mid-infrared wave- 
lengths can break this degeneracy, providing more precise 
measurements of a p lanet's orbital and physical parameters 
ijColon fc Ford! [20091 ) . Observing at these wavelengths also 
allows for reduced differential atmospheric effects, which can 
further improve the quality of the transit light curve. 

The detection of transiting Earth-size planets around 
solar-like stars requires such high photometric precision 
that astronomers had long assumed that characterizing such 
planets could only be done from spa ce to avoid the dele - 
terious effects of Earth's atmosphere (|Borucki et al.l [l985l ^ . 
Narrow-band photometry with large telescopes provides an 
alternative path towards high photometric precision. While 
inefficient for planet searches, near-simultaneous narrow- 
band observations provide additional opportunities for the 
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characterization of known planets (or previously identified 
planet candidates) , such as the measurement of atmospheric 
absorption. 

In this paper we describe a novel observational tech- 
nique for high-precision transit photometry, using near si- 
multaneous observations in multiple narrow bandpasses 
( 92J - We describe our data reduction in JJ]and light curve 
analysis in iJU In ij5]we present results for TrES-2 and TrES-3 
and demonstrate that this technique can achieve a very high 
photometric precision. Finally, in SjHl we compare our results 
to those found in the literature, and we discuss the impli- 
cations for studying the atmospheric composition of giant 
planets and for characterizing the bulk properties of super- 
Earth-size planets, including those in the habitable zone of 
main-sequence stars. 



2 OBSERVATIONS 

The 10.4-m Gran Telescopio Canarias (GTCfl is located 
at the Observatorio del Roque de los Muchachos, on the 
island of La Palma. While the GTC began scientific op- 
erations in March of 2009, commissioning of the telescope 
and first light instruments is ongoing. We describe some 
of the first scientific observations with the GTC using the 
Optical System for Imaging and low Resolution Integrated 
Spectr oscopy (OSIRIS) i n the tunable filter (TF) imaging 
mode (|Cepa et aUbood . l2003h . OSIRIS includes two 2048 
x 4096 pixel E2V 44-82 CCDs which provide a maximum 
unvignetted field of view of 7.8 x 7.8 arcmin. OSIRIS of- 
fers a suite of filters, including a TF which allows the user 
to specify custom bandpasses with a central wavelength of 
651-934.5 nm and a FWHM of 1.2-2.0 nm. Note that the TF 
bandpass is not uniform across the field of view, with the 
effective wavelength decreasing radially outward from the 
centre. Therefore, for the wavelengths used in these obser- 
vations, a difference of ~10 nm exists between the tuned 
wavelength at the optical centre and the wavelength ob- 
served at 4 arcmin from the optical centre (i.e. near one 
edge of the CCD). 

We observed one transit each of two giant extrasolar 
planets, TrES-2b and TrES-3b, using GTC/OSIRIS. The 
field of view was chosen so that: 1) the target and a "pri- 
mary" reference star were observed at the same wavelength 
using the same CCD, and 2) several additional "secondary" 
reference stars were observed with the same CCD but at 
different distances from the optical centre (and thus at dif- 
ferent wavelengths). Due to the restrictions in CCD orien- 
tation, different chips in the OSIRIS CCD mosaic were used 
for the TrES-2 and TrES-3 observations. 

During each transit, observations alternated between 
two bandpasses centred on 790.2 and 794.4 nm at the lo- 
cation of the target. These wavebands were chosen to mini- 
mize atmospheric effects by minimizing water vapor absorp- 
tion and skyglow. Each observation was followed by 33.4-s 
of dead time for readout. We use lxl binning but utilize a 
fast readout mode (500 kHz) in order to decrease the dead 
time between exposures. Recent and future CCD controller 
upgrades will reduce the dead time between exposures. The 
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telescope was also defocused to increase efficiency and to re- 
duce the impact of pixel-to-pixel sensitivity variations. Even 
with a slight defocus, the resulting point-spread functions of 
the stars were fairly well defined (i.e., not doughnut-shaped). 

The TrES-2 (V ~11.4) observations took place under 
photometric conditions on 2009 June 25 (dark time) from 
1:58 to 5:43 UT, during which the airmass ranged from 
~1.07 to 1.44. The defocussed FWHM of the target varied 
from ~1.6-2.8 arcsec (~12.5-22 pixels) during the observa- 
tions, while the actual seeing was ~1.2 arcsec. The auto- 
guider system kept the images aligned within a few pixels, 
with the target's centroid coordinates shifting by less than 
4 pixels in the x-direction and less than 6 pixels in the y- 
direction. We used 80-s exposures, resulting in an overall 
cadence of 3.78 min for each bandpass. Four of the images 
were excluded from our analysis due to tracking problems. 
The last ~16 min of out-of-transit (OOT) data were affected 
by twilight and are not included in the analysis. 

The TrES-3 (V ~12.4) observations took place under 
photometric conditions during grey time, starting at 23:53 
UT on 2009 August 10 and continuing until 2:54 UT the 
next morning. During the observations, the airmass ranged 
from ~L13 to 2.22, the defocussed FWHM of the target 
varied from ~1. 4-2.1 arcsec (~11-16.5 pixels), and the tar- 
get's centroid coordinates shifted by less than 3 pixels in 
the x-direction and less than 5 pixels in the y-direction. The 
exposure time was 120-s, yielding a 5.11 min observing ca- 
dence for each bandpass. Problems with the telescope caused 
our observations to begin shortly after ingress had already 
begun. 



3 DATA REDUCTION 

Standard IRAF procedures for bias subtraction and flat field 
correction were used for the TrES-2 observations. However, 
during preliminary analysis, we found that the total number 
of counts in the dome flat fields (with a given exposure time) 
were decreasing with time. Thus, for flat fielding, we include 
only those dome flats taken after the first 30 min of the 
lamp being turned on, by which time the lamp intensity 
had stabilized. For TrES-2, we combine 40 (out of 60) dome 
flats for each observed wavelength (790.2 and 794.4 nm). 
After summing over either set of flats, the average number of 
counts per pixel was 1.71 xlO 6 . For the TrES-3 observations, 
we perform similar reductions as for TrES-2, but we only use 
45 of the 95 total dome flats taken for each wavelength. After 
summing over multiple flats, the average number of counts 
per pixel was 5.28xl0 5 (1.04xl0 6 ) for 790.2 (794.4) nm. For 
the TrES-3 data reduction, we replace the median bias frame 
with the median of a series of dark frames taken with the 
same exposure time as the observations (albeit taken several 
months after the actual observations took place), since early 
observations revealed a higher than expected dark current. 
We note that subtracting darks was beneficial only for the 
TrES-3 data, based on the RMS scatter of the out-of-transit 
light curve for TrES-3. The opposite is true for TrES-2; i.e., 
the RMS scatter of the out-of-transit light curve was smaller 
when bias frames were subtracted (instead of dark frames). 
This difference is likely due to the fact that the dark frames 
for both targets were taken in October 2009, rather than 
the exact same conditions as the transit observations (i.e., 
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several months after each transit was observed). Therefore, 
each data set was affected differently by its respective set 
of dark frames, resulting in the darks proving to be useful 
only for the TrES-3 observations. Recent hardware upgrades 
have reduced the dark current. 

Due to the TF's small bandpass and position-dependent 
wavelength, all observations contain sky (OH) emission 
rings. Therefore, we performed sky subtraction on all im- 
ages based on the IRAF package TFrecjJ which estimates 
the sky background, including rings due to sky emission. We 
performed aperture photometry on each target and reference 
star using the standard IDL routine APEE0. We repeated 
this analysis using a range of aperture radii, and measured 
the root mean square (RMS) scatter of the flux ratio (target 
over sum of references) outside of transit in each colour. Af- 
ter considering the results in both bandpasses, we adopted 
an aperture radius of 44 pixels (5.6 arcsec) for stars in the 
field of TrES-2 and 38 pixels (4.8 arcsec) for stars in the field 
of TrES-3. Note that the use of TFred removes the need for 
a sky annulus, so we did not input one in the APER routine. 

To account for any atmospheric extinction, we fit lin- 
ear airmass trends to each reference star's light curve (LC; 
computed by dividing the flux for a given star by the sum 
of the flux from all the other reference stars). For each ref- 
erence star, we discarded points that resulted in a flux ratio 
greater than 3<r from the mean (typically one and at most 
two points per reference star). We computed the reference 
flux at each time as the weighted sum of the flux of the 
remaining reference stars. The flux in the target aperture 
was divided by the reference to compute the final LCs. We 
found that using an ensemble of six to eight reference stars 
resulted in a smaller RMS OOT scatter than using just the 
primary reference, despite the fact that the secondary ref- 
erences were observed at different wavelengths and that the 
reference stars in the field of TrES-2 (TrES-3) ranged over 
~2 (3) magnitudes in brightness. Since most of the reference 
stars were fainter than the targets, we conclude that adding 
additional faint references helps to reduce the photon noise 
for a larger ensemble and balance stellar variability within 
different reference stars. 

The uncertainties in the flux ratio were estimated by 
computing the quadrature sum of the photon noise for the 
target and the reference ensemble [median values are 0.331 
(0.339) and 0.244 (0.249) mmag for TrES-2 and 0.426 (0.433) 
and 0.210 (0.214) mmag for TrES-3 for 790.2 (794.4) nm], 
the uncertainty in the sum of sky background and dark cur- 
rent for the target and the reference ensemble [median values 
are 0.0751 (0.0737) and 0.103 (0.105) mmag for TrES-2 and 
0.0998 (0.0840) and 0.0688 (0.0667) mmag for TrES-3 for 
790.2 (794.4) nm], and the scintillation noise for the target 
and primary reference [median values are 0.0627 (0.0630) 
mmag for TrES-2 and 0.0778 (0.0766) mmag for TrES-3 for 
790.2 (794.4) nm]. The uncertainty in the sum of sky back- 
ground and dark current was estimated by performing aper- 
ture photometry on the sky frames produced by TFred at 
the specific location of each target and reference star. The 
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Table 1. Relative Photometry of TrES-2 and TrES-3 



Observed A HJD Relative Flux Uncertainty 



TrES-2 

790.2nm 2455007.5845 0.99963 0.00043 

790.2nm 2455007.5871 0.99969 0.00043 

790.2nm 2455007.5923 1.00006 0.00043 



794. 4nm 


2455007.5858 


0.99956 


0.00044 


794.4nm 


2455007.5910 


0.99996 


0.00044 


794. 4nm 


2455007.5936 


0.99984 


0.00044 




TrES- 


-3 




790.2nm 


2455054.4976 


0.99855 


0.00048 


790.2nm 


2455054.5011 


0.99472 


0.00048 


790.2nm 


2455054.5047 


0.98969 


0.00049 


794. 4nm 


2455054.4993 


0.99800 


0.00049 


794. 4nm 


2455054.5029 


0.99266 


0.00049 


794. 4nm 


2455054.5064 


0.98936 


0.00049 



The time stamps included here are for the times at mid-exposure, 
and the relative flux has been corrected for slopes/trends in the 
light curves (see J5]for more details). 



scintillation n o ise wa s estimated using the rela tion given by 
iDravins et"aD (|l998l l. based on lYound ljl967h . We caution 
that this is an empirical relation and has not been tested 
for large telescopes at excellent sites like La Palma. Thus, 
we consider the expression for scintillation noise to be only 
a rough estimate. Nevertheless, it demonstrates that scin- 
tillation is expected to be only a very small contribution to 
the total error budget, thanks to the narrow filter bandpass. 
Readout and digitization noise as well as flat field noise are 
negligible compared to photon noise and are not included 
in the calculation of the uncertainties in the flux ratio, so 
we include estimat es for these noise sourc es (as based on 
relations given by ISouthworth et all l2009h only for refer- 
ence. The median readout noise is 7.44 electrons per pixel 
or 0.0116 (0.0123) mmag for TrES-2 and 7.35 electrons per 
pixel or 0.0150 (0.0158) mmag for TrES-3 for 790.2 (794.4) 
nm. The median flat field noise is 6.07e-4 (6.07e-4) electrons 
per pixel or 0.00274 (0.00269) mmag for TrES-2 and 1.07e-3 
(7.63e-4) electrons per pixel or 0.00778 (0.00535) mmag for 
Tr ES-3 fo r 790.2 (794.4) nm. Based on the relation given 
bv lHowelll (|2006t ). we find the median total uncertainties in 
the flux ratio for each exposure to be 0.434 (0.443) mmag for 
TrES-2 and 0.495 (0.499) mmag for TrES-3 for 790.2 (794.4) 
nm. We note that our estimated uncertainties are somewhat 
larger than the measured residual RMS given in i]5.2l While 
the estimated and measured precision are marginally con- 
sistent, it is possible that we overestimated the uncertainty 
in individual measurements. 

We perform the above analysis for both bandpasses and 
each target. The resulting photometric time series is re- 
ported and shown (with some corrections; see Sj5} in Table 
[1] (full version available online) and Figures [1] and [2] 
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Figure 1. Normalized light curves (a), residuals (b,c), and colour 
(d) for nearly simultaneous observations at 790.2 and 794.4±2.0 
nm of TrES-2 as observed on UT 2009 June 25. The filled circles 
are observations and the lines in panel (a) show the best-fitting 
models. In panel (a), the 794.4-nm light curve has been arbitrarily 
offset by 0.015. Panels (b) and (c) show residuals from the fits for 
the 790.2 and 794.4 nm light curve. The colour of the residuals is 
shown in panel (d). 



4 LIGHT CURVE ANALYSIS 

Before fitting models to the light curves, we first applied 
the external paramete r deco rrelation (EPD) technique (see, 
e.g., lBakos et al.ll2007l T2010) to each light curve in order to 
remove any systematic trends that are correlated with the 
following parameters: the X and Y centroid coordinates of 
the target on the image frames, the sharpness of the tar- 
get's profile [approximately equal to (2.35/FWHM) 2 ] and 
the airmass. We then performed the following analysis on 
the decorrelated light curves. 

We fit the flux ratio with a standar d planet transit 
mode l that includes a quadratic LD law l|Mandel &: Agoll 
2002). We parametrize the LC model using: the time of 
mid-transit (to), impact parameter (6 = acosi/i?*), tran- 
sit duration (from first to fourth contact; D), planet-star 
radius ratio (p = R p /R+), average OOT flux ratio, a linear 
slope for the OOT flux ratio (a), and two LD coefficients 
(ci = ui + U2 and C2 = ui — W2), where u\ and M2 are 
the lin ear and quadratic LD coefficients of Mandel fc Agoll 
(|2002h . 

We used the publicly available code mpfitfut^ to 
perform Levenberg-Marquardt minimization of \ 2 to iden- 
tify a best-fitting model for the transit photometry. The 
initial guesses for 6, D and p are based on esti mates 
(and their uncertai nties) from iHolman et al.l (120071) and 
ISozzetti et "all J2007|) for TrES-2b and lGibson et all (|2009l ) 
and ISozzetti et all (|2009T ) for TrES-3b. The initial guesses 
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for to were base d on the ephemeris o f iRabus et all (2009) 
for TrES-2 and ISozzetti et all l|2009T i for TrES-3. In each 
case described below, we repeat the local minimization us- 
ing an array of initial guesses for the modeled parameters 
based on the published uncertainties and concluded the re- 
sults of the non-linear model fitting were not sensitive to our 
initial guesses. 

Similarly, we tested several sets of values for ci and 
C2 based on the theoretical LD coefficients estimated for 
each star. Specifically, both LD coefficients were computed 
for our specific bandpasses for a grid of stellar models 
using PHOEBE (|Pr§a fe Zwitterl 120051 ). We interpolated 
in (T e //, log g, [Fe/H]) to estimate ci and C2 for each 
of several sets of s tellar parameters. W e adopted stel- 
lar p arameters from Sozzetti et al.l (|2007T ) for TrES-2 and 
from ISozzetti et~ahl (120091 1 for TrES-3. but we note that 
these abundance measure ments are systematicall y lowe r 
than those determined by lAmmler-von Eiff et"al] (2009). 
IValenti fc Fischerl |2005h emphasize that different studies 
can yield different results due to, e.g., the analysis tech- 
niques and/ or stel lar models used, but bo t h ISozzetti et al.l 
|2007l . 120091 ) and lAmmler-von Eiff et all (|2009h use very 
similar techniques to estimate stellar parameters. Therefore, 
we chose to interpolate LD coefficients over the range of 
the stellar p arameters plus t heir measure ment uncertainties 
as given by ISozzetti et all (20071 , [2009), since their anal- 
ysis used slightly higher res olution spectra than used by 
lAmmler-von Eiff et all (2009). As an extra precaution, we 
confirmed that their measurement uncertainties were real- 
ist ic by comparing them to the adopted uncertainties given 
bv lValenti fc Fischerl l|2005l ). which were based on the anal- 
ysis of observations of the solar spectrum as reflected by the 
asteroid Vesta. 

Most previous studies use two parameter limb darken- 
ing laws and fit for one or both limb darkening parameters. 
ISouthworthl (20081 ) investigated the effects of various limb 
darkening models and assumptions. He concluded that a lin- 
ear limb darkening law was often inadequate, so a non-linear 
law should be used, but that two parameter fits were typi- 
cally highly degenerate. In the cases of TrES-2 and TrES-3, 
two-parameter limb darkening models are particularly de- 
generate, since both have a large impact parameter and the 
planet only probes the s t ellar surface brightness near the 
stellar limb. ISouthworthl (|2008l ) also found that including 
uncertainty in limb darkening parameters was important to 
obtain rea listic uncertainties for transit model parameters. 
Therefore. ISouthworthl (2008) recommended holding C2 fixed 
and fitting for ci . Initially, we tried this approach using sepa- 
rate limb darkening parameters for each star and each band- 
pass and fit for ci . We found that the best-fitting limb dark- 
ening coefficients for one star and our two bandpasses could 
differ significantly, even though stellar atmosphere models 
predict similar limb darkening coefficients in the two nearby 
bandpasses. Therefore, we adopt an alternative approach. 

For both TrES-2 and TrES-3, we model the LCs us- 
ing four different scenarios. First, we fit separate models to 
the flux ratios at the two wavelengths (columns 1 and 2 of 
Table [2]), fixing both ci and C2 at a pair of self-consistent 
values based on PHOEBE models using a single set of spec- 
troscopic parameters. The transit and limb darkening pa- 
rameters in these first two scenarios need not necessarily be 
self-consistent. Comparing these first two models allows us 
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to evaluate the sensitivity of our results to the two different 
bandpasses. Second, we correct each LC based on the previ- 
ous results for the fitted slope and mean OOT flux and fit a 
single model to the the LCs for both bandpasses simultane- 
ously (column 3 of Table |5J. In this scenario, the transit and 
planetary parameters (to, b, D, p) are forced to be the same 
for both bandpasses (self-consistent if we neglect the possi- 
bilities of wavelength-dependent planet radius and contami- 
nation from background light). For this and the subsequent 
scenario, we fix all four limb darkening parameters (two for 
each band pass) at self-consistent values based on PHOEBE 
models. Finally, we fit a single model to the two (corrected) 
LCs, but allow for separate values of D and p for each LC 
(column 4 of Table [2]). For the final scenario, we continue to 
force the fit to have the same to and b for both bandpasses 
since these should be the same regardless of the bandpass 
used for observations, but allow the models in column 4 to 
have different values of p and D in the two bandpasses so as 
to account for a potential difference in the planet radius in 
the two bandpasses or blending of a putative binary star (see 
3H| . Since we have observations of a single transit in mul- 
tiple bands and impose additional constraints that all four 
limb darkening parameters be self-consistent, one would ex- 
pect our method to result in a larger RMS scatter about the 
model than the fitting procedures used by previous authors 
and developed for analyzing transit observations at a single 
waveband. 



The primary differences in the model parameters among 
the four scenarios can b e traced to differe nt sets of limb 
darkening parameters. As lSouthworthl |2008r i demonstrated, 
fixing both LD coefficients at their theoretical values could 
produce measurement uncertainties that are too small. To 
account for the uncertainty in the limb darkening model, we 
repeat each of the analyses varying the spectroscopic param- 
eters by the published uncertainties and thus the calculated 
limb darkening coefficients. We report model parameters for 
the best-fitting model with the smallest \ 2 value, but present 
their parameter uncertainties based on the complete set of 
best-fitting values and uncertainties estimated for all allowed 
LD parameters (see Table [2]). For example, we estimate the 
upper error bar for a given parameter as the result of sub- 
tracting its best-fitting value (as given by the model with 
the smallest \ 2 value) from the maximum sum of a fit value 
and its associated measurement uncertainty, considering the 
full set of models computed for the different LD coefficients. 

To further investigate the measurement e rrors, we ap- 
plied the '"Prayer Bead" method (see, e.g., iDesert et al.l 
2009). Specifically, we construct synthetic light curves by 
calculating the residuals from the initial best-fitting light 
curve model, performing a circular shift, and adding the 
shifted residuals back to the best-fitting light curve model. 
Then, we conducted the above analysis on the synthetic data 
sets. We use the dispersion of the fit parameters to estimate 
the effects of any other systematic noise sources not removed 
by the EPD technique. We discuss the results from this anal- 
ysis as well as the formal 1-a errors (as determined from the 
covariance matrix) for the best-fitting parameters in the fol- 
lowing section for the case presented in column 4. 
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Figure 2. Same as in Figure [T] but for TrES-3 as observed on 
UT 2009 August fO. In panel (a), the 794.4-nm light curve has 
been arbitrarily offset by 0.025. 

5 RESULTS 

We show the flux ratios and the best-fitting models (after 
correcting for slopes/trends based on results given in column 
4 of Table [2} in Figures [T] and [2] Note that the time series 
given in Table[T]was corrected based on results from column 
4 of Table [2] as well. The best-fitting LC parameters for 
TrES-2 and TrES-3 are given in Table [2] and described in 
H5.ll We analyze the residuals and discuss the photometric 
precision in &I5.2I 



5.1 Planetary Parameters 

For the joint analysis presented in column 4 of Table [2] 
(which we adopt as our baseline model), the best-fitting 
parameter values and their formal 1-a uncertainties for 
TrES-2 are: to = 2455007.64605 ± 0.00013 (HJD), b = 
0.8642 ± 0.0036, D (790.2nm) = 0.07365 ± 0.00061 days, 
D (794.4nm) = 0.07481 ± 0.00063 days, p (790. 2nm) = 
0.1253 ± 0.0008, and p (794.4nm) = 0.1245 ± 0.0008. Sim- 
ilarly, for the joint analysis of TrES-3, the best-fitting 
parameter values and their formal 1-a uncertainties are: 
to = 2455054.52447 ± 0.00010 (HJD), b = 0.8356 ± 0.0073, 
D (790. 2nm) = 0.05817 ± 0.00065 days, D (794.4nm) = 
0.05652 ±0.00069 days, p (790. 2nm) = 0.1695 ± 0.0025, and 
p (794.4nm) = 0.1631 ±0.0022. We note that for both TrES- 
2 and TrES-3, these uncertainties are typically smaller than 
those in Table [2l since those also account for uncertainty in 
the LD model. Additionally, the 1-a errors for the median 
best-fitting parameters from the Prayer Bead analysis for 
both TrES-2 and TrES-3 are slightly larger than, but still 
comparable to, the formal 1-a uncertainties, with the largest 
deviations occurring for the uncertainty in the radius ratio. 
We also computed the RMS of the best-fitting parameter 
values from the Prayer Bead analysis as an additional check 
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Table 2. System Parameters of TrES-2 and TrES-3 



Parameter 



Value 



(Ai+-A 2 £ 



(Ai 



A 2 l£ 



TrES-2 



to-2455000 (HJD) 7.6460ll°'o°o2i 
b 8647+ ' 0051 

£»(days) 0.07369±°;°°g*j 

v 1254+ - 0050 

ci (fixed) 0.220 

c 2 (fixed) -0.084 



7 64610+ 00021 
(.O*DlU_ ooi9 

0.8638±°-°°§g 

n n747s+ - 00224 
0.07478_ Q007g 

1245+ - 0054 
u.iz^o_ Q Q010 



0.232 
-0.086 



7 64605 +(U) "" 1:i 
l-"1OVO_ 00013 

8639+ 0036 
u.oojy_ 034g 

07415+ - 00188 
u.U(*J.o_ 00051 

1249+ 0048 
u.±z<U)_ 0007 



Ai : 0.220 
A 2 : 0.232 
Ai : -0.084 
A 2 : -0.086 



-+0.00013 



Ai : 
A 2 : 
Ai 
A 2 



0.8642: 
: 0.07365^ 
: 0.07481 : ' 



7.646051 . 00013 

+0.0030 
0.0349 
+0.00 2 14 
0.00061 
+0.00176 
0.00064 
n l9KQ+0.005 2 
u.izjo_ 000g 

1245+ ' 0046 
u.iz^o_ 000g 

Ai : 0.220 

A 2 : 0.232 
Ai : -0.085 
A 2 : -0.098 



TrES-3 



to-2455000 (HJD) 54.52450io 1 'oooi5 

b 0.8276^5 

^(days) 0.05779 

P 0.1672+°;°™ 

d (fixed) 0.653 

c 2 (fixed) 0.263 



54 52445+ - 00015 

U -«' d < -0.0453 
+0.00 2 2 
0.00086 



0.05554 



1641+ ' 0067 

U.1D^1_ 00g7 



0.222 
-0.086 



+0.0001 2 



54.52449_o.oooio 

n 871 fi+ ' 0123 
0.8716_ 0455 

-0.00160 

-0.00061 



0.05608 



1662+ ' 0046 

U.1DDZ_ Q04g 



Ai : 0.220 
A 2 : 0.222 
Ai : -0.081 
A 2 : -0.086 



54 52447+ - 00011 

04-0^41 -0.00011 
~+0.048 2 
-0.0073 



0.8356 1 



Ai : 
A 2 : 
Ai 
A 2 



+0.0007 2 



0.05817_ 00203 

05652+ 00069 
u - uoooz -0. 00174 

■ 1695 +0 ' (,(,4G 

• u - ±Dyo -0.0049 

■ 1631+ 0048 

• u - 1DJ1 -0.0044 



Ai 


0.653 


A 2 


0.625 


Ai 


0.263 


A 2 


0.233 



a We label the bandpass centred on 790.2-nm as Ai for simplicity. 
b We label the bandpass centred on 794.4-nm as A 2 for simplicity. 

c The best-fitting parameters were determined from the joint analysis of the two light curves for each target, 
with only the LD coefficients determined separately for the two wavelengths. See text for additional details. 
^ Same as in column 3, but also with D and p determined separately for each wavelength. 



on our measurement errors. We find that accounting for the 
uncertainty in the LD model sufficiently accounts for the 
distribution of the best-fitting parameters as derived from 
the Prayer Bead analysis. 

For both TrES-2 and TrES-3, the best-fitting values for 
each parameter are consistent between the different models 
presented in Tabled but only after accounting for uncertain- 
ties in the limb darkening model. The most notable differ- 
ence in our results with different modeling procedures is for 
TrES-3 when comparing the impact parameter in columns 
2 and 3 (analysis for the 794.4-nm light curve and the joint 
analysis with common planet parameters) with the other 
two models for TrES-3. This is at least partly due to the 
near degeneracy between impact parameter and limb dark- 
ening model. Indeed, the LD coefficients for the models in 
column 2 and 3 are significantly different from the best val- 
ues for the other two cases. While we were able to place 
tight constraints on the impact parameter for TrES-2b, the 
precision for TrES-3b is significantly reduced, since we are 
modeling an incomplete light curve and have no data prior 
to ingress. 

The best-fitting slopes and their formal uncertain- 
ties as measured in the individual LCs before perform- 
ing the joint analyses for the TrES-2 and TrES-3 LCs are 
0.00012 ± 0.00172 and -0.00020 ± 0.00180 days -1 for the 
790.2 and 794.4 nm LCs for TrES-2, and 0.00256 ± 0.00501 



and 0.00092 ± 0.00540 days -1 for TrES-3. While the esti- 
mated slopes for TrES-2 are fairly consistent, the difference 
in the estimated slopes is much larger for TrES-3 than for 
TrES-2. Note that the measured slopes account for the ef- 
fects of both differential atmospheric extinction (due to some 
reference stars being observed at a different wavelength than 
the target and primary reference) and real astrophysical 
variability in the colour of the target and/or reference stars. 
Our observations are not generally able to disentangle the 
two potential causes of a slope. If the slope were primarily 
due to differential extinction, then observations of multi- 
ple transits should produce consistent results. On the other 
hand, if stellar variability is significant, then each transit 
would need a separate slope parameter. In the case of TrES- 
3, the large difference in slopes between the two bandpasses 
suggests that differential extinction is not responsible. How- 
ever, we caution that the measurement of the slope around 
transit is likely affected by the lack of pre-transit data. Thus, 
we emphasize the importance of acquiring complete transit 
LCs as well as extended, uninterrupted OOT baseline data. 

For TrES-2, the uncertainties in the OOT flux (as de- 
termined from the individual LC analyses) were ~ 0.000091 
(0.000097) for the 790.2 (794.4) nm observations, based on 
26 (25) exposures totaling ~35 (33) minutes of integration 
time. The uncertainties for the OOT flux for TrES-3 are 
slightly larger, perhaps due to the lack of pre-transit obser- 
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vations and covariance with the slope. However, the TrES-3 
observations were more stable overall and allowed for a sin- 
gle baseline of OOT data that was longer than either the 
pre- or post-transit baselines for the TrES-2 observations. 
Both observations demonstrate the capability for very high- 
precision measurements that could enable the detection of 
super-Earth-sized planets around solar-like stars and/or the 
characterization of the atmospheres of giant planets orbiting 
bright stars (see $Qlfel6)l. 



5.2 Light Curve Residuals 

We computed the residual flux by subtracting the best- 
fitting models given by the joint analysis presented in col- 
umn 4 of Table [2] from the data (see Figures [Tp,c and [2)3, c). 
Note that some of the residuals show evidence of additional 
systematic noise sources that were not removed with the 
EPD technique. While we do not know the origin of these 
systematics, we note that the uncertainty estimates from the 
Prayer Bead analysis are consistent with our best-fitting pa- 
rameters (as discussed above, in 35.1) 1. 

We estimate the photometric precision to be 343±45 
and 412±43 parts per million (ppm) for observations of 
TrES-2 at 790.2 and 794.4 nm, based on the root mean 
square deviation of the residuals. Despite high airmass (up 
to ~2.22), the TrES-3 observations produced a precision of 
470±64 and 424±57 ppm at 790.2 and 794.4 nm. Because 
the estimated precisions for both targets are smaller than (or 
comparable to, based on the upper limits computed from the 
expected standard deviations) the estimated measurement 
uncertainties (see Sj3j) , we take them to be consistent with 
the theoretical limit. 

Due to the long, uninterrupted OOT baseline obtained 
for TrES-3, we take the TrES-3 observations as the best- 
case scenario to achieve high precisions consistent with the 
theoretical limit and present the standard deviation of the 
time-binned residuals for the TrES-3 LCs in Figure [3] The 
residuals for the 790.2-nm LC are consistent with the trend 
expected for white Gaussian noise, while the residuals for the 
794.4-nm LC deviate at binning factors larger than 5. When 
binning the uncorrelated 790.2-nm residuals over longer ob- 
serving times (~40 min), we estimate a precision of 146 ppm, 
which is sufficient to detect the transit of a super-Earth-size 
planet (Js 1.3i?®) orbiting a solar-like star or an Earth-size 
planet orbiting a star smaller than 0.7 Rq (assuming ~120 
min of observations between 2nd and 3rd contact and 3-rr 
confidence level). 



5.3 Transit Colour 

In Figures [TJi and[2ji, we present the colour of the residual 
fluxes. The colour was computed by averaging each pair of 
points in the 790.2-nm LC residuals and dividing by the 
flux residuals of the 794.4-nm LC. We find no significant 
colour deviation in either of the TrES-2 or TrES-3 systems. 
We estimate RMS precisions of 519 ppm and 502 ppm for 
TrES-2 and TrES-3, which are both slightly larger, but still 
consistent, with the precisions estimated for the individual 
LCs. We found no significant difference in the colours in and 
out of transit for TrES-2. For TrES-3, there was a slight slope 
in the colour during transit, but we do not consider this to be 



0.001 on 



0.0001 




Binning Factor 



Figure 3. Standard deviation of time-binned residuals for TrES- 
3 as a function of the number of data points per bin (n). The 
blue square (red triangle) symbols are the standard deviations 
of the 790.2 (794.4) nm binned light curve residuals. The solid 
curves show the trend expected for white Gaussian noise 
Each exposure was separated by 5.11 minutes, thus the present 
observations only allow us to test timescales of up to ~40 minutes 
due to the limited duration of the observations. 



significant given the large uncertainty in the slopes fit to the 
individual LCs due to the lack of pre-transit observations. 



6 DISCUSSION 

The first exoplanet observations from the GTC provided 
excellent photometric precision, despite problems with the 
telescope, high airmass and relatively poor atmospheric con- 
ditions. 

Our results for TrES-3 are c onsistent wit h pre- 
vious high-precisio n observations l| Gibson et al.l 120091 ; 
ISozzetti et al.l 12009), Our measurements of planet-star ra- 
dius ratio, impact parameter and transit duration bracket 
those of previous studies. Our measured transit time has 
an uncert ainty of ^9.5- s and occurs 44-s after that pre- 
dicted by [Sozzetti et al.l 1120091) and 14-s after the up dated 
ephemeris of I Gibson et all (|2009l 'l. lGibson et~all l|2009h con- 
ducted a transit timing analysis of TrES-3b and they find no 
evidence of transit timing variations. Given the differences 
between various ephemerides and the possibility for stellar 
variability to contribute additional timing noise, we do not 
consider such a difference significant. 

Our results for for Tr ES-2 b are also in good agr eement 
with iHolman et all (|2007f l and ISozzetti et all (|2007n . How- 
ever, we note that all reported values for the radius ratio 
could be affected by a binary companion to TrES-2, which 
has a separation of 1.089 ± 0.008 arcsec fr om the primary 
star a nd A i' ~ 3.7mag as reported by iDaemgen et all 
(2009). While the companion star is clearly included within 
our aperture, our results do not account for the flux of the 
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companion in our analysis. In this case, our model param- 
eter p serves as a depth parameter which can be related 
to the actual pla net-star radius ratio for a given amount 
of contamination |Daemgen et al.| [2009). If we assume that 
the primary-secondary star flux ratio in our bandpasses is 
similar to that in i' , then the planet-star radius ratio would 
increase by ~1.6% relative to our estimates in Table [2] In 
principle, the different colour of a blended star could result 
in different transit depths when observed in multiple differ- 
ent bandpasses. While our results are suggestive of such a 
transit depth difference, the difference is not statistically sig- 
nificant. Given the high-precision of our observations and a 
possible feature in the in-transit data (around mid-transit), 
we suspect that the apparent differences in p may be due 
to variations in the stellar surface brightness^ In light of 
these potential complications, near-simultaneous photome- 
try in multiple bandpasses could be particularly useful for 
recognizing when a companion or background object is af- 
fecting the photometry. 

Several recent studies have considered the possibility of 
vari ations in the transit times and/or durations of TrES- 
2b dMislis fc Schmittl 120091; iMislis et al.l l2010l : iRabus et ail 
l2009l : lRaetz et al.l2009l ; IScuderi et a.l.ll2O10h 7 Our best-fitting 
transit time has a measurement uncer tainty of ^11-s and i s 
offset from the predicted ephemeris of lHolman et al.l (120071') 



i prec 

bv_~4 minutes, IRabus et al.l d2009h by_~142-s, iRaetz et all 
(|2009h by ~16-s and IScuderi et ail I^Old ) by ~64-s. Given 
the differences in various ephemerides, observing bandpasses 
and limb darkening models, we do not consider the tran- 
sit time offset to be significant evidence for transit timing 
variations. This is consi s tent w ith the r ecent transit timin 
analysis of IRaetz et ail (|2009h . While IRabus et all ([200' 
raised the possibility of sinusoidal variations in the transit 
ephemeris due to an exom oon, they do not find re sults of 
the magnitude proposed bv lMislis fc Schmittl (I2009T I. 

Due to TrES-2b's large impact parameter, the transit 
duration is very sensitive to changes induced by additiona l 
planets or other bodies i n the system |Miralda-Esc udc 2002). 
IMislis fc Schmittl (|2009l ) and IMislis et all (|2010l ) claim that 
the transit duration decreased by three minutes between 
2006 and 2008, and they argue that a third body is a natu- 
ral explanation for this change. If we include uncertainty due 
to unknown limb darkening models, we c annot definitively 
rule ou t th e shorter durations s uggested bv lMislis fc Schmittl 
( 2009) and IMislis et all |20ld ). When interpreting putative 
differences in transit duration/impact parameter, one should 
be mindful of dif ferences in the vario us filters and the stellar 
LP m odel used |Scuderi et al.ll201ol ). Indeed. IScuderi et alj 
(|2010h report no significant change in the orbital inclina- 
tion/transit d uration when they compa re their observations 
with those of lO 'Donovan et alj (|2006h (taken in the same 
filter). Since our observations use unique bandpasses, it is 
not (yet) possible to compare them to previous observations 
with the same bandpass. Since TrES-2 is in the Kepler field, 
Kepler observations should soon shed light on this matter, 



5 Indeed, previous studies have found other indications of pos- 
sible stellar activity, i n cluding variations in the OOT flux 
llO'Donovan et al.l 1 20061 ; IRaetz et all |2009|). Our observa tions 
show no signs of a putative "second dip" jRaetz et al]|2009h . but 
only extend for ~1 hour after the end of transit. 



as the observations will cover longer timescales and therefore 
multiple transits all within the same passband. 

In summary, we find that large ground-based obser- 
vatories are capable of achieving high-precision differential 
photometry. Systematic effects due to variable atmospheric 
extinction can be minimized by the use of a small far- 
red or NIR bandpass chosen to avoid sky absorption lines. 
Thus, large, ground-based observatories can help character- 
ize super-Earth-size planets discovered by ongoing transit 
searches, including planets in the habitable zone of main- 
sequence stars for which transits last several hours. We 
anticipate that this innovative technique for high-precision 
photometry will enhance the ability of current and future 
large optical/NIR observatories to study the properties of 
both Earth-like planets (e.g., size, density and orbit) and 
giant planets. For example, the high photometric precision 
could allow for the detection of the occultation of hot transit- 
ing giant planets like TrES-3b at optical wavelengths. While 
the occultation depth probes the temperature and dynamics 
of the planet's atmosphere, the time and duration of occul- 
tation provide a powerful probe of the orbit's eccentricity. 

The technique is even more powerful when observa- 
tions in multiple bandpasses can be obtained nearly simul- 
taneously and thus during the same transit. Such obser- 
vations could help verify that the transit of a candidate 
planet is not due to stellar variability plus rotation by ob- 
serving in multiple wavelengths that reduce the effects of 
stellar variability and/or enhance the contrast of spots or 
plage. Similarly, nearly simultaneous observations in mul- 
tiple narrow bandpasses could be used to characterize the 
atmospheres of giant planets. For example, HST and HET 
have used transmission spectroscopy of HD 209458b and HD 
189733b to measure Na I absorption in the planetary at - 
mospheres j|Charbonneau et al.ll2002l ; iRedfield et all 120081 ') . 
Future GTC observations could be used to perform simi- 
lar measurements of atmospheric features such as Na I and 
K I absorption in planets transiting bright stars, especially 
once planned improvements in the OSIRIS CCD controller 
software reduce the dead time between exposures. 

Considering the increasing numbers of known transit- 
ing exoplanets, we hope that this novel technique will be 
an additional useful tool to help characterize a variety of 
planetary properties. 
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